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The most recent evidence for geodetic precession observed in binary radio pulsars 
is presented and discussed. It is demonstrated how an analysis of these results 
can be used to study theories of gravity, stellar evolution and pulsar emission 
theory. In order to highlight the observational strategies, an overview of the applied 
techniques is given. 

1 Introduction 

The discovery of pulsars in 19670 marked the beginning of a new era in astrophysics. 
Besides studying pulsars as fascinating astrophysical objects in their own right, 
these sources serve as invaluable tools in the research of a large variety of physical 
and astrophysical problems. Studied subjects range from plasma and solid-state 
physics under extreme conditions to investigations of the interstellar medium. One 
can study stellar and binary evolution as well as high precision astrometry and 
planetary ephemerides. Applications are found in cosmology, quantum physics and, 
of course, in gravitational physics. The latter application became possible with the 
discovery of the first binary pulsar by Hulse & Taylor in 19740. As they immediately 
realized, a pulsar in a binary orbit represents a highly stable and accurate clock 
orbiting in the gravitational field of (usually) compact stars. They were fortunate 
to discover in PSR B1913+16 the first double neutron star system. These systems 
are doubly rare because they are formed less often and also since they are more 
difficult to detect (significant Doppler smearing of the pulse signal occurs when the 
orbits are less than about 1 day). 

Meanwhile many more binary pulsars have been discovered. At the time of 
writing about 5%, or 70 pulsars out of more than 1300 known, are in binary systems. 
Three pulsars are in orbit with normal or giant stars, and two pulsars are members 
of planetary systems. The rest of the binary pulsar population has compact stars, 
i.e. white dwarfs or neutron stars, as companions. The majority are orbited by 
white dwarfs of typically 0.2Afq,-| .although there is a small but growing population of 
heavy white-dwarf companions&crEI. Pulsar-white-dwarf systems have been used to 
perform unique tests of certain aspects of gravitational theories^, for instance testing 
the existence of gravitational dipole radiation as predicted by many n alternative 
theories of gravityu or the violation of the strong equivalence principlecl. However, 
for tests of theories of gravity in a strong field limit, observations of double neutron 
star systems still yield the most exciting results. 

Only five double neutron star systems (DNSs) are known. Orbital periods range 
from 7.7 hours for PSR B.1913+160 to 18.8 days in the case of the most recent dis- 
covery, PSR J1811 — 17360. Presently, two of the five DNSs allow the determination 
of more than two post-keplerian (PK) parameters. The PK parameters describe 
relativistic corrections to Newtonian physics, represented by the five standard Ke- 
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plerian parameters, and can be written as functions of the pulsar and companion 
masses, m p and m c , and the Keplerian parameters. These functions will differ for 
different theories of gravity, but when Keplerian and PK parameters are measured 
by observations in a manner independent of a particular theory of gravity, the ob- 
servations can be compared to theoryEj. If two PK parameters are measured, a 
given theory of gravity will produce values for the pulsar and companion mass. 
More generally, a measurement of n PK parameters describes n curves in a two di- 
mensional m p -m c plane whose shape and position depend on the theory of gravity 
being applied. If the chosen theory is an accurate and adequate description of the 
physics, all n curves should meet in a single point. A measurement of at least three 
PK parameters hence allows one to perform tests of gravitational theories 

The measurement of PK parameters has been possible thanks to a powerful 
technique called pulsar timing. When giant radio telescopes are used to measure 
the arrival time of photons emitted by the pulsar, a breath-taking accuracy can be 
achieved in experiments which are impossible to perform on Earth or in the solar 
system. For details, we refer to the contribution, by Stairs (these proceedings) and 
to reviews like that of Backer & Hellings (1986)E2l. Despite the remarkable successes 
achieved with pulsar timing, one particular aspect of gravitational physics can only 
be studied adequately if the structure of the received pulses itself is analysed care- 
fully. 

In general relativity, the proper reference frame of a freely falling object suffers 
a precession with respect to a distant observer, called geodetic precession. A direct 
measurement of this effect is the scientific goal of the Gravity Probe-B satellite ex- 
periment (see Everitt, these proceedings). In a binary pulsar system this geodetic 
precession leads to a relativistic spin-orbit coupling, analogous to spin-orbit cou- 
pling in atomic physics. As a consequence, the pulsar spin precesses about the 
total angular momentum, changing the relative orientation of the pulsar towards 
Earth. In such a case, we should expect a change mite radio emission received from 
the pulsar, as first proposed by Damour & Ruffinita very soon after the discovery 
of PSR B19I3+16. Hence, in contrast to "only" accurately registering the arrival 
times of pulsar signals, we are even more interested in the emission properties when 
we want to study geodetic precession. 

In this paper we will summarize how to study geodetic precession in binary 
pulsars. As we have to analyse the emission properties of radio pulsars, we first 
provide a brief introduction into pulsars and their radiation properties. We report 
observations for PSR B1913+1 6 and discuss models used to describe detected pulse 
shape changes. We demonstrate how information about the system geometry can 
be used to learn about the effects of asymmetric supernova explosions and "kicks" 
imparted to newly born neutron stars. Observations of two more binary pulsars, 
PSRs B1534+I2 and Jl 141— 6545 are discussed as they provide more evidence for 
geodetic precession. 



2 Pulsars 
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Table 1. Selected parameters for the pulsars discussed. Uncertainties in the last digits are quoted 
in brackets. 





PSR B1913+16 1 


PSR B1534+12 2 


PSR J1141-6545 il 


P (ms) 


59.029997929613(7) 


37.9044403237164332(10) 


393.897833900(3) 


Pb (d) 


0.322997462736(7) 


0.4207373013796296135(10) 


0.197650966(6) 


e 


0.6171308(4) 


0.2736777(5) 


0.171881(9) 


a sini (lt-s) 


2.3417592(19) 


3.729463(3) 


1.859470(14) 


m p (M ) 


1.4411(7) 


1.339(3) 


? 


m c (M Q ) 


1.3879(7) 


1.339(3) 


? 



Baylor & Weisberg (1989), 2 Stairs et al. (1998), 3 Kaspi et al. (2000) 



2.1 The Clock 

Pulsars are rotating neutron stars. From timing measurements of binary radio 
pulsars we determine the masses of pulsars to be within a quite narrow range of 
1.35 ± 0.04MqEJ. Model calculations involving different equations of state produce 
results for the size of a neutron star quite similar to the very first calculations 
by Oppenheimer & Volkoff (1939)E-3, i.e. about 20 km in diameter. Such sizes 
are consistent with independent estimates derived from modelling light-curves and 
luminosities of pulsars observed in X-rays, e.g. for PSR J0437— 4715E3. 

These massive and compact stars rotate with typical periods of about 0.6 s for 
"normal" pulsars, and periods as low as 1.56 ms for millisecond ("recycled") pulsars 
(cf. Fig. [l(]). Considering that this smallest rotational period for PSR B1937+21 
corresponds to more than 38000(!) revolutions per minute makes the large stability 
of pulsars as clocks easy to understand. 

Typical magnetic field strengths at neutron star surfaces are of order 10 12 G. 
By comparison, the Earth's magnetic field is only 1 G. We can infer these field 
strengths from the slow increase in pulse period due to the loss of rotational energy 
by the emission of magnetic dipole radiation, or more directly from synchrotron 
lines in X-ray spectra of X-ray pulsarstS. 

2.2 The Clockwork 

The tick of the pulsar clock is provided by a narrow radio beam centred on the 
magnetic axis of the pulsar, which is inclined to the rotation axis (see Fig. |). If 
the pulsar beam is directed towards Earth once per rotation, radio telescopes on 
Earth may be able to pick up a pulsed radio signal. This periodic beacon sent 
by the pulsar clock is usually rather weak due to the distance of the pulsar, but 
also due to the small size of the emission region, which we expect to be only a few 
hundred km acrosstll. Nevertheless, the radio emission is actually rather intense as 
we measure brightness temperatures exceeding 10 30 K, which we can only explain 
by assuming a coherent emission mechanism. After 35 years of pulsar research, the 
details of the actual emission process still elude us, but at least we have some basic 
understanding sufficient to perform the experiments described later. 

As the neutron star rotates with its magnetic field, an induced electric 
quadrupole field pulls out easily charges from the stellar surface (the electrical force 
exceeds the gravitational force by a factor of 10 12 !), surrounding the pulsar with 
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Figure 1. left) A pulsar is a rotating, highly magnetised neutron star. A radio beam centred on 
the magnetic axis is created at some distance to the pulsar. The tilt between the rotation and 
magnetic axes makes the pulsar in effect a cosmic lighthouse when the beam sweeps around in 
space; right) while individual pulses vary in shapes and strength (zoomed in, top), average profiles 
are stable (full profile, bottom). The typical pulse duty cycle is only ~4%. 



dense plasma. The magnetic field forces the plasma to co-rotate with the pulsar 
like a rigid body. This co-rotating magnetosphere can only extend up to a distance 
where the co-rotation velocity reaches the speed of lightQ. This distance defines the 
so-called light cylinder which separates the magnetic field lines into two distinct 
groups, i.e. open and closed field lines. Closed field lines are those which close 
within the light cylinder, while open field lines would close outside. The plasma 
on the closed field lines is trapped and will co-rotate with the pulsar forever. In 
contrast, plasma on the open field lines can reach highly relativistic velocities and 
can leave the magnetosphere, creating the observed radio beam at a distance of a 



"Strictly speaking, the Alfven velocity will determine the co-rotational properties of the magne- 
tosphere. 
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few tens to hundreds of km above the pulsar surface (see Fig. [l]). 
2.3 The ticks 

Individual Pulses and Average Pulse Profiles 

Individual pulses reflect the instantaneous plasma processes in the pulsar magneto- 
sphere at the moment when the beam is directed towards Earth. The dynamics of 
these processes results in often seemingly random individual pulses, in particular 
when viewed with high time resolution (see Fig. [I]). Despite this variety displayed 
by the single pulses, the mean pulse shape computed by averaging a few hun- 
dreds to few thousands of pulses is incredibly stableta. In contrast to the snapshot 
provided by the individual pulses, the average pulse shape, or pulse profile, can 
be considered as a long-exposure picture, revealing the global circumstances in the 
magnetosphere. These are mostly determined by geometrical factors and the strong 
magnetic field, leading to very stable pulse profiles. Apart from a distinct evolution 
with radio frequency, the same profiles are obtained, no matter where and when 
the pulses used to compute the average have been observed. [] 

Beam Shapes 

The observed pulse profiles show a large variety of shapes. Although each pulsar 
exhibits a slightly different profile - almost like a unique fingerprint - a systematic 
pattern can be recognized. The most simple .model successfully describing the 
beam shapes is that of a hollow cone of emissiorcj. It is based on the idea that the 
outermost open field lines, which show the largest curvature among the "emitting 
field lines" , should be associated with the strongest emission, leading naturally 
to a cone-like structure. Observations show that this picture is oversimplified, 
since we also often observe components near the centre of the cone, a so called 
core component, in particular at lower frequencies. However, the combination of a 
hollow cone with a core compouerit can be successfully applied to understand the 
vast majority of pulsar profilesEJEaEj. Depending on the way how our line-of-sight 
cuts the emission cone, different pulse profiles are observed (Fig. ||). 

Polarisation — Signatures of Geometry 

An important and most useful property of pulsar radio emission is its typical high 
degree of polarisation. The radiation is often 100% elliptically polarized, and it 
is usual practice to separate the polarisation into linearly and circularly polarized 
components. The linear component is often the far dominating one, although pul- 
sars with circular components as strong as 30% or more are not uncommon. The 
degree of polarisation decreases with increasine-irequency, leaving pulsars more or 
less completely de-polarized at high frequencies^. At low frequencies, however, the 
polarisation serves as a useful diagnostic tool, both to measure the magnetic field 

b For completeness we should note that there is a very small fraction of pulsarsQIEJB exhibiting 
profile changes on short time scales of minutes or hours. In these cases, the profiles do not exhibit 
random shapes either, but the pulsar seems to switch between a few, usually two, distinct profiles. 
The origins of these phenomena are not understood. On secular time scales, these pulse profiles 
are still stable. 



kramer: submitted to World Scientific on February 1, 2008 




Figure 2. left) Hollow cone model to explain the observed pulse profiles as different cuts through 
the emission beam; right) very often an additional "core" component is observed near the centre 
of the beam, in particular at low frequencies. 



of the ionized interstellar medium and, in particular for our purposes, to obtain 
information about the viewing geometry. 

The strong coupling of the outwards-moving plasma to the magnetic field lines 
in the pulsar magnetosphere has the effect that the plane of polarisation of the 
linear component is determined by the plane embedding the corresponding field 
line. The observed position angle (PA) of the linear polarisation is then given by 
the projection of this direction onto our line-of-sight. The result is an S-like curve 
of the PA whose shape depends on the angle between the rotation and magnetic 
axes as well as on the distance of our line-of-sight to the magnetic pole. If our line- 
of-sight cuts the emission beam close to the magnetic axis, the PA changes rapidly 
when crossing the centre of the beam. If we are cutting the cone further away from 
the pole, the transition is much smoother and the PA swing much flatter. 

By measuring the polarisation characteristics of a pulsar, we can in principle 
win information about tha pulsar's orientation towards us. In practice, fitting 
this rotating vector modeSl (RVM) often turns out to be difficult. Although the 
majority of observed PA swings can be well described by the RVM after correcting 
for sometimes occurring orthogonal modes (i.e. jumps of the PA by nearly 90° which 
are probably magnetospheric propagation effects), the uncertainties in the obtained 
angles representing the geometry are typically large. The reason is not a failure 
of the model, but simply the small size of beam of most pulsars, which provides 
constraints to the fit for only the small fraction of the pulse period when the pulse 
is actually observed, which is typically only 4%. 

The geometrical origin of the PA swing across the pulse is impressively demon- 
strated by its independence upon observing frequency]^. Even if the results of a 
RVM fit to the PA data are not very well constrained, an observed change in the 
swing curve will immediately indicate a change in viewing geometry. 



c after correcting for propagation effects due to the interstellar medium, i.e. Faraday rotation 
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Figure 3. left) Pulsar emission is often highly polarized. The PA of the linear component usually 
performs an S-like swing across the profile, right) the change in PA can be understood as a 
projection of the field line direction onto our line-of-sight. 



3 Evidence for Geodetic Precession — Part I. 



Immediately after the discovery of the 59-ms binary pulsar PSR B1913+16 it was 
realized that the system should exhibit a measurable amount of geodetic precession 
if the p^aar-spip axis is misaligned with respect to the orbital angular momentum 
vector E3H'ESEj We will discuss the possible origin for a misalignment of the spin 
vectors later and first concentrate on the expected observational consequences. 



3.1 Precession rate 



The precession rate is given byl 




O - f—\ rji2/3 m c (Am p + 3m c ) 1 

Up -{pJ ' 1q ' 2(m p +m c )^ ' l-e* UJ 

where Pb is the period and e the eccentricity of the orbit. We express the masses m p 
and m c in units of solar masses (Mq) and the define the constant T Q = GMq/c 3 = 
4.925490947/xs. G denotes the Newtonian constant of gravity and c the speed of 
light. 

The values for the orbital period and the eccentricity of the PSR B1913+16 
system can be obtained from timing observations (see Table |]) . These observations 
have also resulted in the measurements of three PK parameters, providing the 
most stringent test of gravitational theories in the strong field limit so farl3E3E-3. 
General relativity has passed all these tests with flying colours. We therefore have 
safe and accurate measurements for the masses of both pulsar and its companion: 
vrip = 1.4411 ± 0.0007 and m c = 1.3879 ± 0.0007. (Strictly speaking these observed 
masses differ from the intrinsic ones by an unmeasurable Doppler factor due a radial 
velocity, which cannot be determined from timing measurements. We can safely 
neglect this small difference here.) Using these values we obtain a precession rate 
of f2 p = 1.21 deg yr _1 . Since the orbital angular momentum is much larger than 
the pulsar spin, the orbital spin practically represents a fixed direction in space, 
defined by the orbital plane of the binary system. Given the calculated precession 
rate, it takes 297.5 years for the pulsar spin vector to precess around it. 
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Figure 4. Precession changes the cut of our line-of-sight (indicated by circle) through the emission 
cone. The consequence should be a change in the width and possibly also shape of the pulse profile. 
Since the distance between line-of-sight and magnetic axis is changing, an alteration of the position 
angle swing is expected. 



3.2 First Signs - First Riddles 

As a result of the precession the angle between the pulsar spin axis and our line- 
of-sight should change with time, so that different portions of the emission beam 
are observed. Consequently, one expects changes in the measured pulse shape, in 
particular in the profile width, as a function of time (see Fig. |J). In the extreme 
case, the precession may move the beam out of our line-of-sight and the pulsar may 
disappear from the sky until it becomes visible again. 

Since the precession rate predicted by general relativity is reasonably large, the 
pulse profile was naturally studied closely in opier to detect possible changes. Fi- 
nally, Weisberg et al. (1989, hereafter WRT89)E3, reported a change in the relative 
amplitude of the two prominent profile components (Fig. |^, left). Analyzing ob- 
servations made at 1408 MHz from 1981 to 1987, they determined a change in the 
amplitude ratio of 1.2 ± 0.2% yr _1 , i.e. a weakening of the leading component. 

These first signs of profile changes could be indeed explained by a changing cut 
through the emission beam and hence as the first evidence for geodetic precession. 
However, if the emission beam exhibits an overall hollow-cone like shape, one would 
also expect a change in the separation of the two components rather than only a 
change in relative intensity. This was not observed. WRT89 determined an upper 
limit of AW < 0.°06 for a change in component separation in the given time 
interval. In order to reconcile this constant profile width with the seen amplitude 
changes, they argued that the beam structure is irregular ancLpatchy rather than 
cone-like, i.e. similar to models by Lyne & Manchester (1988)E3. 

With different cuts through the emission beam, the distance of our line-of-sight 
to the magnetic axis should also change with time. A change in the PA swiug 
would be expected. Cordes, Wasserman & Blaskiewicz (1990, hereafter CWB90)E3 
therefore studied polarisation data of PSR B1913+16 to compare profiles and PA 
swings obtained at frequencies between 1397 and 1416 MHz from 1985 to 1988. 
CWB90 neither detected very clear changes in the pulse shaped], nor could they 



^Presumably, this result was caused by using profiles measured at different frequencies, allowing 
frequency evolution to confuse the measurements. 
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Figure 5. left) Pulse profile of PSR B1913+16 as presented by Taylor & Weisberg (1989). The 
change in component amplitudes studied by Weisberg & Taylor (1989) is clearly visible, right) 
profile of B1913+16 in 1995 and component amplitude ratio as presented by Kramer (1998). 

find any significant change in the PA swing. CWB90 pointed out that their results 
and those of WRT89 could be still consistent with a hollow cone emission beam. 
They proposed that the existence of a core component, which is prominent at lower 
frequencies!^ (see Fig. ||), is causing the change in the component amplitude ratio. 
Finally, they suggested that the lack of detected changes in profile width and PA 
swing was due to a special precession phase at the time of observation. 

The situation remained somewhat unsatisfactory. While the amplitude changes 
detected by WRT89 seemed to be clear indications of the presence of geodetic 
precession, the actually expected behaviour, i.e. a change in profile width and PA 
swing, was not observed. 

3.3 The Times, They are a Changing 

In 1994 the pulsar B1913+16 was included in a timing program at the 100-m radio 
telescope at Effelsbcrg, which was initiated to take over the regular timing of mil- 
lisecond pulsars during the upgrade of the 300-m Arecibo telescope in Puerto RicpEa. 
The analysis of profiles measured between 1994 and 1998 by Kramer (1998)c3 re- 
vealed that the profile components were still changing their relative amplitude, 
consistent with the rate first determined by WRT89 (see Fig. ||). Even more in- 
teresting, however, was the final detection of a small but significant change in the 
separation of the components (see Fig. |6|). In order to model the long-expected 
decreasing width of the profile, two simple assumptions were made, i.e. those of a 
circular hollow cone-like beam and a precession rate as predicted by general rela- 
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tivity. Both well justified assumptions (see above) lead to a model which has only 
four free parameters: the misalignment angle A between the pulsar spin and the 
orbital angular momentum, the inclination angle between the pulsar spin axis and 
its magnetic axis, a, the radius of the emission beam, p, and the precession phase 
given by the reference epoch Tq. With these parameters the component separation 
at an epoch t is given bytiil 



W{t) = 4 sin" 1 



/sin 2 (p/2) - sin 2 ((Jft) - cp/2) 
smas'm8(t) 



(2) 



The angle 5 measures the distance of our line-of-sight to the pulsar spin axis at 
the closest approach to the magnetic pole. Due to the precession, this distance will 
change with time: 

cos 5(t) = cos A cos i + sin A sin i cos <f>{t) (3) 

where <p(t) is the precession phase, 

4>{t) = il p ■ (To - t) (4) 

With the PK parameters measured by pulsar timing, general relativity allows one 
to compute the value of sin?, i.e. the sine of the orbital inclination angle. In a 
fortuitous edge-on geometry of the pulsar orbit, we would even be able to measure 
this value as the shape of a Shapiro delay visible in the timing data (cf. PSR 
B1534+12 below). For PSR B1913+16, we simply compute a value of i = 47.°2, or 
equivalently i = 180 — 47.2 = 132. °8 (since we only obtain sinz from timing). 

The best fit of this model to the data allows four equivalent solutions. One pair 
of solutions corresponds to i = 47. °2, the other pair to i = 132. °8, respectively. 
The remaining choice is given by the unknown relative orientation of the pulsar 
spin and the orbital angular momentum, i.e. as to whether the pulsar rotation is 
pro-grade or retro-grade. As we will see later, it can be argued that a retro-grade 
case is less likely. This leads to the following best solutionEJ 

i = 132.°8, a = 153tf°, A = 22±f°, p = 9.0tt° and T = 2128 ± 4 

as this is the only . one .which gives the right sense of PA swing when compared with 
polarization data.E3e3 In fact, this system geometry, which is solely determined 
from the_change in component separation, is fully consistent with the results of a 
RVM fitc3. Moreover, the obtained-misalignment angle of A = 22^g is in excellent 
agreement with earlier simulationsc3 made to study the number of observable DNSs 
which predicted A « 20° as a typical value for PSR B1913+16-like systems. 

The obtained best fit shown in Fig. || leads to the prediction that the pul- 
sar will disappear from the sky around the year 2025! Moreover, it also implies 
that the component separation remains almost unchanged for about 60 yr, which 
corresponds to a likelihood of 20% to observe the pulsar in that phase given the 
precession period of about 300 yr. Looking back, it is now easy to understand why 
WRT89 were not lucky to detect changes in the component separation. Similarly, 
computing the change in PA swing which had to be measured by CWB90 for a 
positive detection of a geometry change, produces a value which is only slightly 
larger than their estimated detection limit. 
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Figure 6. left) Coordinate system used to describe the performed modeling by Kramer (1998); 
right) measurements of the component separation as a function of time including the best fit 
(taken from Kramer 1998). 

About ten years after the first indications of geodetic precession, it is finally 
possible to provide solid evidence for its existence by modelling profile changes. 
Most exciting, however, may be the prediction of the pulsar's disappearance in 
about 25 years. Shortly before this, the leading component wilL disappear if it 
continues to weaken with the measured rate (cf. Istomin 1991)e-l Reappearing 
again around the year 2220, PSR B1913+16 will, in total, only be observable for 
about a third of the precession period (cf. Figs. ^| & ||). 

3.4 Further Results and Update 

After the completion of the Arecibo upgrade, Weisberg & Taylor (2000)0 obtained 
new measurements, confirming the results reported above. Thanks to the superior 
sensitivity of the Arecibo telescope, they not only measured the same decrease in 
component separation, but could also measure a general decrease in profile width 
at several intensity levels (Fig. |7]l_.The derived misalignment angle of A = 14(±2)° 
with an upper bound of A = 22°c3 is in good agreement with the fit to the Effels- 
berg data. The obtained data quality even allows one to do the following exciting 
experiment: Since our line-of-sight moves through the emission beam, each profiles 
represents a slightly different cut through the beam structure. By adding measured 
profiles in the right order, it is possible to reconstruct a real 2-D map of a pulsar 
beam for the first time (Fig. [?]). The first results of this mapping procress are 
surprising: Given the profiles available so far, the pulsar beam seems to be elon- 
gated. Although the possibility of a non-spherical pulsar beam has been discussed 
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Figure 7. left) Decrease of PSR B1913+16's profile width at different intensity levels as determined 
by Weisberg & Taylor (2000) using the Arecibo telescope, right) resulting 2-D pulsar beam map 
reconstructed from the different profiles. 



several times in the literature, even with claims for elongation in both latitudinal 
and longitudinal directions, convincing evidence had not been presented. 

During their data reduction, Weisberg & Taylor (2000) perform a brilliant step 
to separate the measured profiles into odd and even parts. While the even profiles 
are symmetric to a chosen midpoint, the odd profiles contain all the asymmetries. 
The addition of both odd and even parts then reproduces the actual pulse shape. 
This trick removes the complicating effects of the change in relative component 
amplitude from the symmetric profiles, so that these can be used to map the 2-D 
beam structure by applying a specific mapping function. 

Since an elongation of the pulsar beam would have interesting implications (for 
instance for birth rate calculations), it is stimulating to compare the results of 
Weisberg & Taylor (2000) with the outcome of an independent study that employs 
reverse engineering (Kramer 200f , in prep.). Starting with a circular cone, cuts 
through the emission beam are computed for different lines of sight. Using the re- 
sulting profiles, a beam map is computed, following the same procedure as applied 
by Weisberg & Taylor (2000). We are interested to see whether we are able to 
reproduce the observed profiles without resorting to elongated beams. Detailed re- 
sults of this study will be presented elsewhere, but one interesting effect can already 
be obtained: in order to produce the asymmetric profiles as observed in reality, a 
core component is placed in the beam at a slightly off-centre position. The core is 
visible at low, frequencies (see Fig. |^), and its possible influence was already pointed 
out earliercZlEZl. Using the profiles simulated for given epochs, we can determine 
both the expected component amplitude ratio and the component separation and 
compare them with the observations. A snapshot of these simulations is shown in 
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Figure 8. Further evolution of PSR B1913+16's profile in a simple hollow cone model including 
an off-centre core component which is visible in the 430-MHz profile (left, Taylor & Weisbcrg 
1982). The right figure is a snapshot of a simulation measuring component amplitude ratio and 
separation. While the component amplitude ratio first decreases as observed, we can predict that 
it will rise again as soon as the line-of-sight leaves the core dominated region. 



Fig. pi As a result, we can not only reproduce the observed values, but we can even 
make a prediction: when the line-of-sight leaves the central region dominated by 
the core, the amplitude ratio should increase again and may eventually reach unity 
again (Kramer 2001, in prep.). 

In Fig. |^ we present the latest measurements obtained with the Effelsberg tele- 
scope. We also add Weisberg & Taylor's previously unpublished Arecibo data. As 
visible from the new measurements, the previously observed trend continues. In the 
coming years, the decrease in profile width should accelerate drastically, eventually 
allowing us to turn this qualitative test of general relativity into a quantitative one. 
For demonstration, we attempted to treat the rate of precession as a free parameter 
rather than keeping it fixed to the value predicted by general relativity. The result, 
ftp = 1.2 ± 0.2 deg yr" 1 , (shown in Fig. ^|) is by no means a precision test and 
not even a proof of principle! In order to succeed, one has to be confident that the 
beam structure is modelled correctly. Certainly, this will be the most difficult part 
of the job. However, it may be possible that an iterative process will be finally 
successful. Moreover, we did not make much use of the available polarisation in- 
formation yet. While we used PA data to exclude certain classes of models, we did 
not invoke_RVM fits because of the mentioned large uncertainties. Karastergiou et 
al. (2000)l3 presented a method where polarisation and pulse shape data can be 
combined without relying on RVM fits which will become very useful. 

Combining the polarisation data will be a crucial test, as we should definitely 
detect changes in the PA swing finally. Eventually, one may obtain reliable measure- 
ments of fi p . However, the most reliable test would be, of course, the re-appearance 
of the pulsar at about 240 years after it disappears — if we can wait that long. 
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Figure 9. Latest results for measurements of the still decreasing component amplitude ratio (left) 
and the decreasing component separation (right). The left plot contains unpublished data from 
Jodrell Bank, while the right plots includes previously unpublished data by Weisberg & Taylor 
(2000). During the fit leading to the results shown on the right, we treated Q p as a free parameter. 



4 Asymmetric Supernova Explosions 

Geodetic precession occurs if the pulsar spin is misaligned with the orbital momen- 
tum vector. For PSR B1913+16, we can measure this misalignment angle using 
the observed precession effects to obtain A « 22°. The question about the origin of 
this angle is far from trivial. In fact, a non-zero A is the imprint of the violent birth 
event which created the DNS. It is the ultimate proof that supernova explosions 
are asymmetric. 

4.1 The Birth of a DNS 

In the beginning, a DNS starts off as a massive binary system where the more mas- 
sive companion eventually explodes and becomes a pulsar. The pulsar is probably 
born with an initial period of 10-30 ms and a high spin-down rate P (upper left 
region in P — P-diagram, Fig. As the rotation slows down to a period of a 
few seconds, the pulsar emission eventually ceases after tens of million years or so, 
moving the pulsar to the lower right corner of the P — P-diagram. When the mas- 
sive companion evolves and overflows its Roche-Lobe, the binary system reaches 
a High-Mass X-Ray Binary phase. The dead pulsar accretes matter and angular 
momentum. This spin-up "recycles" the pulsar, which re-appears as a millisecond 
pulsar in the lower left part of the P — P-diagram. The spin period achieved during 
the recycling process depends on the duration of the mass transfer. For low-mass 
systems this phase can be rather long, resulting in periods of a few ms as usually 
observed for pulsars with white-dwarf companions. In the case of high-mass systems 
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Figure 10. left) P — P-diagram as known at the time of writing. Binary pulsars are marked with 
circles. Lines of constant age and surface magnetic field strengths are indicated; right) evolutionary 
stages of the PSR B1913+16 (a) before and (c) after the second supernova explosion. Before the 
explosion all spin vectors are aligned. After the explosion, imparting a kick to the companion, the 
new orbit is misaligned to the previous orbit. 

however, this phase is ultimately limited by the event of the supernova explosion 
(SN) of the companion. 

Just prior to the explosion of the companion, the binary system is expected 
to consist of a recycled pulsar in a circular orbit about its companion, which has 
evolved to a He-staiH. All spin vectors are aligned as a result of the angular- 
moment transfer during the accretion process and a following common envelope 
phase. When the He-star explodes, the pulsar itself is not directly affected by the 
explosion, since its cross section is too small. In contrast, the survival of the binary 
system depends on a number of lucky circumstances. 

If the explosion is symmetric, the survival will be mainly determined by the 
total mass loss. Reality, however, is more complicated, since there is convincing 
evidence that SNe are in fact asymmetric. In that case, a momentum "kick" is 
imparted to the newly born neutron star. Unless_the magnitude and direction of 
the kick are favorable, the binary system disruptsEJ. 

4-2 The Kick — Lessons from Geodetic Precession 

The actual process causing the kick is still poorly understood. Most of the pro- 
posed models, which can be roughly classified as electromagnetic kicks, neutrino 
driven kicks or hydrodynamical kicks, either produce velocities too small compared 
to observations, or require unrealistic conditions such as extremely large magnetic 
fieldsEil. It turns out that it is rather difficult to produce kick velocities in excess 
of, say, 300 km s _1 . In a recent work, Wex et al. (2000) have used the system in- 
formation obtained by geodetic precession for PSR B1913+16 to derive constraints 
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Figure 11. Profiles of PSR B1534+12 as measured by Stairs et al. (2000) at 430 MHz (left) and 
1400 MHz (right). Changes due to geodetic precession are clearly visible. 



on the kick mechanis: 




Since the recycled pulsar is unaffected by the explosion of its companion, its 
spin vector will continue to point in the direction of the pre-SN orbital angular mo- 
mentum. The new orbital plane, and hence the new orbital angular momentum, is 
determined by the direction and magnitude of the kick and will be usually different 
from the old orbit. As a result, the two vectors are misaligned and gravitationally 
coupled - geodetic precession occurs(Fig. |l0| ). 

It is possible to write down a set of equations relating the pre- and post-SN 
orbit using conservation of momenta, conditions for the orbit, typical sizes of He- 
star etcEl Knowing that all spin vectors before the SN had been aligned, and 
assuming that the system velocity after the second SN is dominated by the effects 
of the kick, we can infer those kick parameters that are necessary to explain the 
post-SN system characteristics. This procedure is possible for PSR B1913+16 (Wex 
et al. 2000)Ej: Because of its nature as a DNS, both stars do not interact, but the 
system's evolution is solely determined by the emission of gravitational waves and 
its motion in the galactic potential. Both can be calculated backwards in time. The 
(justified) assumption that the system was born in the galactic plane, marks the 
termination point of this backwards computation. The derived orbit and system 
velocity must then be similar to the conditions just after the SN. 

With the unmeasurable radial velocity as the only free parameter, only three 
different cases are consistent with the age of the pulsar and the system's current 
motion away from the plane. Considering all these cases and possible pro- and 
retro-grade solutions, Wex et al. (2000) reach a number of interesting conclusions 
among which are: a) on the average, a larger kick is required to explain a retro-grade 
solution; b) kick velocities between 300 and 1100 km s^ 1 are necessary; c) indepen- 
dent of radial velocity and considered cases, the kick imparted to the companion 
must have been almost perpendicular to the orbital spin. Given the alignment of all 
pre-SN spins, the kick was also perpendicular to the rotation axis of the exploding 
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star. That has very interesting implications. Firstly, the duration of the kick must 
have been short compared to the rotational period of the exploding star (otherwise 
this velocity component would have disappeared). Secondly, if we assume that the 
magnetic field of the companion was more or less aligned with the rotational axis, 
as it is the case for the Earth or Sun, the kick was also perpendicular to the mag- 
netic axis. The latter rules out all kick mechanisms which have a preference in the 
direction of the magnetic field, as for instance some neutrino oscillation modelsEi 
Yet again, pulsars can obviously serve as diagnostic tools for studying a rather 
different branch of astrophysics. By carefully investigating this particular binary 
system, we are able to derive important constraints for any kick mechanism pro- 
posed, i.e. every proposed model should be able to produce a short kick perpendic- 
ular to the rotation axis. 

5 Evidence for Geodetic Precession — Part II. 

5.1 Geodetic Precession and PSR B 1534+ 12 

The next obvious source to look for evidence of geodetic precession is the sec- 
ond DNS discovered in the disk of our Galaxy, PSR B1534+12E3. This 38-ms 
pulsar happens to be brighter than PSR B1913+16, exhibiting also a much nar- 
rower pulse peak which both lead to a much better timing accuracy. The fortunate 
edge-on orientation of the orbit allows the measurement of two PK parameters de- 
scribing the Shapiro delay of the pulse arrival times. In total five PK parameters 
can be measuredE-3 (see Stairs, these proceedings). With the mass estimates of 
m p = 1.334 ± 0.002 and m c = 1.344 ± 0.002 (Table Q), the expected precession 
rate of n p = 0.51° yr^ 1 is much-Lower than that of PSR-B1913+16 (Fig. [[§). 
Nevertheless, Arzoumanian (1995)E3 and Stairs et al. (2000)E3 succeeded in detect- 
ing secular changes in the pulse profiles providing further evidence for geodetic 
precession acting in binary pulsars (Fig. pi] ). 

The profiles of recycled pulsars with periods, less then 50 ms do not seem to fit 
the classical hollow-cone beam model anymoreo. This may stem from the extreme 
compactness of millisecond pulsar magnetospheres, whose sizes scale with the pulse 
period. The difficulties in modelling the profile changes observed for PSR B1534+12 
in a manner done for PSR B1913+13 are hence unfortunate but not too surprising. 
The situation can possibly be saved by the excellent polarisation data available. 
Since emission can be seen for almost 80% of the pulse period, RVM fits provide 
information about the viewing geometry of PSR B1534+12 with an accuracy much 
better than for most other pulsarsEj. 

5.2 PSR J II4I-6545 — Further Evidence? 

The latest relativistic binary discovered in the Parkes Multibeam Survey!, PSR 
J1141— 6545, was first considered to be a DNS. The short orbital period of 4.7 
hours and an orbital eccentricity of e = 0.171 (Tab. [I]) fit well into the range known 
for DNSs. Soon after its discovery, the first PK parameter could be measured, 
i.e. the periastron advance, u>. Assuming that its value is completely determined by 
relativistic effects, general relativity then allows one to determine the total mass of 
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Figure 12. left) Precession rates as predicted by general relativity for the sources discussed here 
and a further DNS PSR B2127+11C. Unfortunately, its pulsar is rather weak. The precession 
rate for PSR J1141— 6545 is given as a function of the unknown orbital inclination angle; right) 
profiles of PSR J1141-6545 at 430 and 1400 MHz (Kaspi et al. 2000). 



the systemtj, M — m p + m c = 2.4. Assuming a canonical pulsar mass of m p ~ 1.4 
leaves only m c w 1 for the companion. These estimatesJead to the conclusion that 
the companion is almost certainly a heavy white dwarfi 

The exact mass distribution in the system can be determined if a second PK 
parameter is measured, which should be easily possible in the future. For the 
moment, we can make use of the mass function of the binary system which indicates 
that the inclination angle must be close to i = 90° to produce a pulsar mass of 
uip = 1.4. Using these values, the predicted precession rate is the largest for any 
known system, i.e. fl p w 1.35 deg yr . For smaller (though unlikely) inclinations, 
the value even increases (see Fig. |l^ ). 

Obviously, even as a pulsar-white dwarf system, PSR J1141— 6545 is a prime 
candidate to exhibit geodetic precession. It turns out that the same spot of sky had 
been surveyed at least twice ten years ago, but no pulsar was discovered, neither 
at 430 nor at 1400 MHz. Today, the pulsar is fairly strong at both frequencies, 
with a flux density several times that of the detection limit of the previous surveys 
(Fig. |l2|). Unless other effects like man-made radio interference are responsible for 
the previous non-detections, it is conceivable that the pulsar is indeed precessing 
and has moved its beam into our line-of-sight during the past ten years. If that 
is the case, future observations should show profile changes soon. Since the pulsar 
is strong and also polarised, we can hopefully present further detailed evidence for 
geodetic precession in the near future. 

6 Outlook 

It is clear from the presented observational evidence that general relativity correctly 
predicts the existence of geodetic precession. While the presently available data 
provide only qualitative tests so far, quantitative tests may be possible if the pulsar 
emission beam can be modelled sufficiently. Alternatively, polarisation data may 
produce reliable geometry information via RVM fits, so that a measurement of Q p 
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will be possible. 

Although the analysis of pulse structure data is the most approachable method 
to investige geodetic precession, recently Doroshenko et al. (in prep.) investigated 
the effects on pulsar timing in more detail. It turns out that geodetic precession 
may produce after all a measurable effect in the timing. Applying their calculations 
to PSR B1913+16 leads however to the conclusion that this effect may be absorbed 
in the fits for P and P, since the pulsar may only be visible for 60 to 80 years. It 
has yet to be determined as to whether other pulsars could be useful. 

As we have seen, one major consequence of geodetic precession can be the (tem- 
porary) disappearance of a known pulsar from the sky. Over a time span that is 
long compared to typical precession periods, i.e. several hundreds to thousands of 
years, more DNSs will be discovered and measured than without geodetic preces- 
sion. However, for the typical life time of a scientist, the relevant numbers will not 
change: the number of binary pulsars that will disappear from the sky before we 
can detect them should be similar to the number of sources turning up for the first 
time. In other words, we live in a "steady-state precessing pulsar universe" , so that 
birth rate and hence detection calculations for gravitational wave detectors are un- 
affected by the discussed phenomena. While this consideration certainly applies in 
general, it is a very reasonable idea to occasionally re-search the position of known 
DNSs with large precession rates to check for the appearance of previously unseen 
companions. You never know... 
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